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Abstract 

We analyze two different confining mechanisms for Levy flights in the pres- 
ence of external potentials. One of them is due to a conservative force in the 
corresponding Langevin equation. Another is implemented by Levy-Schrodinger 
semigroups which induce so-called topological Levy processes (Levy flights with 
locally modified jump rates in the master equation). Given a stationary probabil- 
ity function (pdf) associated with the Langevin-based fractional Fokker-Planck 
equation, we demonstrate that generically there exists a topological Levy process 
with the very same invariant pdf and in the reverse. 

PACS numbers: 05.40. Jc, 02.50.Ey, 05.20.-y, 05.10. Gg 

1 Motivation 

Levy flights stand for a nickname of a class of symmetric stable stochastic processes. 
These non-Gaussian jump-type processes are not yet deeply accommodated within the 
general statistical physics conceptual and technical imagery. One obstacle is rooted in 
purely technical reasons - a shortage of explicit analytic solutions, e.g. probability den- 
sity functions and transition probability densities. Other obstacles may be attributed 
to somewhat exotic features of Levy flights, like an existence of arbitrarily small jumps 
or that of fat tails of the probability density which preclude the existence of moments 
(in the least of the second one). 

There is also a physically motivated obstacle: no physical thermalization mech- 
anisms have ever been proposed for Levy flights. Moreover, their physical "reason" 
(origin of noise) appears to be exterior to the physical system, with no reliable kinetic 
theory background, and therefore no fluctuation-dissipation response theory could have 
been been set for any stable noise. 

To the contrary, the noise "reason" is definitely an intrinsic feature of the environment- 
particle coupling in case of the standard Brownian motion, based on the kinetic theory 
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derivations. All traditional fluctuation-dissipation relationships find their place in the 
Brownian framework. None of them has been reproduced in the context of Levy flights. 

Nonetheless, Levy flights remain an active research area, with a definite tendency 
to expand their field of applicability. A fraction of papers, relevant to our problem 
(taming Levy flights via suitable confining potentials) is reproduced in the reference 
list, P-P3J. 

In the present paper we address an apparent discord between two families of jump- 
type processes which can be be associated with the same type of the free (Levy noise. 
Namely, most of current research is devoted to Langevin equation-based derivations, 
where a deterministic force is perturbed by Levy (white) noise of interest, [3]- [13]. How- 
ever, in a number of publications, another class of jump-type processes was introduced 
under the name of topologically induced super-diffusions, [H]-[IH]- As we indicate in 
below they belong to a broader class of processes governed by Levy-Schrodinger semi- 
groups. In a different context, they were introduced and investigated in some detail in 
Refs. E], c.f. also 0. 

Fractional Fokker-Planck equations, inferred respectively for Langevin and topo- 
logical processes, appear not to describe the same dynamical pattern of behavior. The 
statement of [16] was that topologically induced fractional Fokker-Planck equation does 
not portray a situation equivalent to any of standard (e.g. Langevin-based) fractional 
Fokker-Planck equations and thus does not correspond to any Langevin equation. 

Confined Levy flights have received some attention in the recent literature [8]-|13j. 
On the other hand, in the discussion of topological Levy processes main emphasis has 
been put on their super-diffusive behavior with some neglect of confining effects and 
the resultant emergence of non-Gibbsian stationary probability densities, [T4]-[T8]. 

We address the latter issue and set general confinement criterions for analytically 
tractable of Cauchy noise-driven processes. To this end, some ideas have been adopted 
from the general theory of diffusion-type stochastic processes where an asymptotic 
approach towards equilibrium (stationary density) is one of major topics of interest, 
[22]. 

To handle topological Levy processes we invoke a theoretical framework of the 
Schrodinger boundary data problem (whose validity is well established in the Brownian 
context, [HI EI])- We shall demonstrate that topologically induced processes of [H]- 
[T8] form a subclass of its solutions with a properly tailored dynamical semigroup and 
its (Feynman-Kac) potential, [SI E]- That allows to take advantage of the existing 
mathematical theory of Levy processes and Levy-Schrodinger semigroups, [3], 0] and 
[SI El [7] , where free Levy noise generators are additively perturbed by suitable confining 
potentials. The theory works well for both Gaussian and non-Gaussian processes. (For 
the record we indicate that, in the Brownian case, the Schrodinger problem incorporates 
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the well known transformation of a Fokker-Planck equation into a generalized diffusion 
equation, [23"].) 

In Section 2 we give a necessary background for the forthcoming analysis of random 
motions in terms of dynamical (Schrodinger and next Levy-Schrodinger) semigroups. 

In Section 3 we extend the Schrodinger boundary data problem Section 2 to Levy 
flights. In Section 4 we address in detail a number of model cases related to the Cauchy 
flights, with emphasis on their confinement mechanisms. 

2 Brownian inspiration: Schrodinger semigroups 

2.1 Gibbsian asymptotics of Smoluchowski processes 

We begin from a brief resume of standard Smoluchowski processes and their general- 
izations. Since in the present discussion we assume forward drifts to be of the gradient 
form, our considerations will be limited (albeit with no substantial loss of generality) 
to one space dimension. 

Let us consider x = b(x,t) + A{t) with (A(s)) = 0, (A(s)A(s')) = V2D5(s - s'). 
Here, b(x, t) is a forward drift of the process, admitted to be time-dependent, unless we 
ultimately pass to Smoluchowski diffusion processes where b(x, t) = b(x) for all times 
and x G R. 

Given an initial probability density po{x). We know that the diffusion process drives 
this density in accordance with the Fokker-Planck equation d t p = DAp — V (b ■ p) . 
We introduce an osmotic velocity field u = DV In p, together with the current velocity 
field v = b — u. The latter obeys the continuity equation dtp = — Vj, where j = v ■ p 
has a standard interpretation of a probability current. 

Presently we pass to time-independent drifts of the diffusion process, that are in- 
duced by external (conservative, Newtonian) force fields / = — VV. One arrives at 
Smoluchowski diffusion processes by setting 

h= —«=-^R VV - W 

mp mp 

This expression accounts for the fully-fledged phase-space derivation of the spatial 
process, in the regime of large f3. It is taken for granted that the fluctuation-dissipation 
balance gives rise to the standard form D = kBT/mfl of the diffusion coefficient. 

Let us consider a stationary asymptotic regime, where j — > = 0. We denote an 
(a priori assumed to exist, [22]), invariant density p* = p*(x). Since i>* = 0, there holds 

&* = = DVlnp* . (2) 

However b = f/m/3 by its very definition does not not depend functionally on the 



probability density. Thus b = b* and therefore p*(x) = (l/Z) exp[—V(x)/kBT]. Our 
outcome has the Gibbs-Boltzmann form. 

Denoting F* = — fc^TlnZ, where Z = J exp(— V/kBT) dx, we have 

p*0) = exp ([F* - y(x)]/A; B T) = exp[2$(x)] . (3) 

Here, to comply with the notation of [191 E] and the forthcoming discussion of a topo- 
logical generalization of the Brownian motion and next Levy flights [H]-[l8], we have 

1/2 

introduced a new potential function $ such that pj = exp $ and b = 

2.2 Schrodinger semigroups and their integral kernels 

Given Eq. (3), let us re-define a probability density of the Smoluchowski process as 
follows 

p(x, t) = 9 m (x, t) exp[$(ar)] . (4) 

This allows to map a Fokker-Planck equation into an associated generalized heat equa- 
tion [23]. Indeed, if p{x,t) obeys the Fokker-Planck equation whose drift b(x) has the 
gradient form, the inferred 8*(x,t) obeys the parabolic equation 

dA = DA9* - V6* , (5) 

provided a compatibility condition 

V(x) = I + Vb] (6) 



2 \2D 
holds true. 

We have here arrived at a contractive semigroup dynamics = [exp(—tH)8 *](0) , 
with the dynamical (Schrodinger) semigroup operator exp(-tH) whose generator —H 
involves the familiar (here, l/2mD rescaled) Schrodinger Hamiltonian 

H = -DA + V. (7) 

The initial t = data read 9*(x, 0). 

For V = V(x),x G R, typically expected to be continuous and bounded from below, 
the operator H = —DA + V is essentially self-adjoint on a natural dense subset of L 2 , 
and the integral kernel (s < t) 



k(y, s, x, t) = (exp[-(t - s)H\j (y, x) 



of the dynamical semigroup exp(-tH) is strictly positive. Assumptions concerning an 
admissible potential may be relaxed. The necessary demands are that H is self-adjoint 
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and bounded from below. Once we have an integral kernel, the semigroup dynamics 
takes the form 9*(x, t) — J #*(?/, s) k(y, s, x, t). 

The key role, in the transformation of the Fokker-Planck operator and equation into 
Schrodinger-type operator and (parabolic) equation, is played by the potential V. We 
recall that —DA is a semigroup generator for the free Brownian motion which yields 
the heat equation d t p — DAp. 

Brownian noise perturbations of conservative force fields may be encoded on two 
equivalent ways. One approach directly employs the Langevin equation, c.f our discus- 
sion of Smoluchowski processes. It is the deterministic force field ~ —W which one 
additively perturbs by white noise. 

Another approach is based on the classic Feynman-Kac formula which uniquely 
assigns a dynamical semigroup kernel (and dynamics) to a concrete potential V(x): 



the conditional Wiener measure i.e. along all sample paths (Brownian trajectories) of 
the Wiener process X(t) which connect y with x in time t — s; X(t) stands for the 
random variable of the Wiener process. 

2.3 Schrodinger 's boundary data problem 

Schrodinger semigroups naturally appear in the framework of so-called Schrodinger 
boundary data problem, [19] . We assume to have a priori given two probability densities 
p(x, 0) and p(x, T), T > and attempt to construct a Markov process that interpolates 
between these boundary data. 

The major outside input is a definite choice of a bounded strictly positive and 
continuous integral kernel k(x, s,y,t), for all times in the interval (possibly infinite) 
< s < t < T. Our considerations are restricted to Feynman-Kac kernels and so to 
Schrodinger semigroups. Each semigroup is selected by prescribing the Feynman-Kac 
potential V, c.f. Eqs. (7) and (9), which needs to be a continuous function. For each 
concrete choice of the semigroup kernel, cf . [J9J [20] . 

By denoting (notice that we have here combined both forward and backward prop- 
agation scenarios in the time interval [0,T]): 




(9) 



In the above we integrate the exp[— J V(X(u), u)du] weighting factor with respect to 




(10) 
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where functions f(x) and g(x) are suitable inital/terminal data, it follows that for all 
< s < t < T there holds 



p(x, t) = 9(x, t)9*(x, t) = J p(y, s, x, t)p(y, s)dy . (11) 

Here 

9(x,t) 

p{y, s, x, t) = k(y, s, x, t) ' (12) 

0{y,s) 

is the transition probability density of the pertinent Markov process that interpolates 
between p(x, 0) and p(x,T). We note that 9*(x, 0) = f(x) while 9(x,T) = g(x). 

In case of the free evolution, by setting 9(x,t) = 9 = const, we effectively trans- 
form an integral kernel k of the ^(R) norm-preserving semigroup into a transition 
probability density p of the Markov stochastic process, e.g. the familiar heat kernel. 

Our 9 and 9* are defined up to constant factors. This freedom is limited by the 
demand that 9-9* = p actually is a probability density of the pertinent Markov process. 

We emphasize that in addition to the semigroup dynamics, to secure the unique- 
ness of an interpolating Markov process, we need an auxiliary motion rule, c.f. [T9] . 
(propagating backwards in time the prescribed terminal function 9(x,T)) 

d t 9 = -DA9 + V9 . (13) 

For time-independent 9 = 9(x) = exp$(x) we note that V = DA9/9. Taking into 
account the definition of the forward drift b = 2DV9/9, [19], we end up with the 
previous compatibility condition (6). 

In case of Smoluchowski processes we encounter asymptotic invariant densities, [22] . 

1/2 

Accordingly, 9 ~ pj and thus 

V = fl%. (14) 

P* 

This result derives from the semigroup argument alone and does not rely on the Gibb- 
sian form of p*. 



3 Levy-Schrodinger semigroups 
3.1 Stable noise and its generator 

The Schrodinger boundary data problem is amenable to an immediate generalization 
to infinitely divisible probability laws which induce contractive semigroups (and their 
kernels) for general Gaussian and non-Gaussian noise models, [5] and [6]-[T5]. A sub- 
class of stable probability laws contains a subset that is associated in the literature 
with the concept of Levy flights. 



6 



Instead of the semigroup generators proper, we shall employ the rescaled Hamilto- 
nians (true generators have an opposite sign) of the form H = F(p), where p = — iV 
stands for the momentum operator (up to the scaled away h or 2mD factor), and for 
— oo < k < +00, the function F = F(k) is real valued, bounded from below and locally 
integrable. Then, for a function f(x) in the domain of the semigroup operator, we have 

[exp(-tH)f](x) = [exp(-tF(p))f(p)] v (x) (15) 

where the superscript V denotes the inverse Fourier transform, and / stands for the 
Fourier transform of /. Let us set 

k t = 

The action of exp(-tH) can be given in terms of a convolution (i.e. by means of an 
integral kernel k t = k(x — y,t) — k(y, 0, x, t)): 

exp(-tH)f = [exp(-tF(p))f(p)] v = f * k t (17) 

where 

(f*g)(x):= [ g(x-z)f(z)dz. (18) 

J R 

We shall restrict considerations only to those F(p) which give rise to positivity 
preserving semigroups: if F(p) satisfies the celebrated Levy-Khintchine formula, then 
kt is a positive measure for all t > 0. The most general case refers to a combined 
contribution from three types of processes: deterministic, Gaussian, and the jump- 
type process. 

We recall that a characteristic function of a random variable X completely deter- 
mines a probability distribution of that variable. If this distribution admits a density we 
can write E[exp(ipX)] = f R p(x) exp(ipx)dx which, for infinitely divisible probability 
laws, gives rise to: 

F(p) = - [ + \exp{ipy) - 1 - -^]u(dy) (19) 

where v{dy) stands for the appropriate Levy measure. The corresponding non-Gaussian 
Markov process is characterized by 

E[exp(tpX t )} = exp[-*F(p)] . (20) 

Accordingly, the contractive semigroup generator follows (keep the minus sign in mem- 
ory) from: F{p) = H. 



[exp{-tF(p)X 



(16) 
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For the sake of clarity we restrict further considerations to non-Gaussian random 
variables whose probability densities are centered and symmetric, e.g. a subclass of 
stable distributions characterized by 

F(p) = \\ p \i>=>H = A|A|^ /2 . (21) 

Here p < 2 and A > stands for the intensity parameter of the Levy process. To 
account for the interval < \i < 2 boundaries, one should rather employ (— A) M//2 
instead of |A| M / 2 , since —A is a positive operator. 

The fractional Hamiltonian H, which is a pseudo-differential operator, by con- 
struction is positive and self-adjoint on a properly tailored domain. A sufficient and 
necessary condition for both these properties to hold true is that the probability density 
of the Levy process is symmetric, [3j. 

The associated jump-type dynamics is interpreted in terms of Levy flights. In 
particular 

F(p) = X\p\ -> H = F(p) = A|V| = A(-A) 1/2 (22) 

refers to the Cauchy process, see e.g. [HI El Ej. 

The pseudo-differential Fokker-Planck equation, which corresponds to the fractional 
Hamiltonian (22) and the fractional semigroup exp = exp(-A|A|^ 2 ), reads 

d tP = -X\A\^ 2 p, (23) 

to be compared with the heat equation d t p = +DAp. 

3.2 Response to external potentials: stationary densities 
3.2.1 Langevin scenario 

In case of jump-type (Levy) processes a response to external perturbations by conserva- 
tive force fields appears to be particularly interesting. On the one hand, one encounters 
a widely accepted reasoning, [8]- [13] where the Langevin equation, with additive de- 
terministic and Levy "white noise" terms, is found to imply a fractional Fokker-Planck 
equation, whose form faithfully mimics the Brownian version, e.g. (c.f. [5], see also [7]) 

x = b(x) + A"(t) => d t p = -V(6 • p) - \\Af 2 p (24) 

where we keep the notation b = f /m/3, f = — W of Eq. (1). 

We emphasize a difference in sign in the second term, if compared with Eq. (4) of 
Ref. [8]. There, the minus sign is absorbed in the adopted definition of the (Riesz) 
fractional derivative. Apart from the formal resemblance of operator symbols, we do 
not directly employ fractional derivatives in our formalism. 



The validity of (24) and temporal details of an approach towards an asymptotic 
invariant density were investigated for the Cauchy-Ornstein-Uhlenbeck proces [7]. That 
safely extends to more general stable (symmetric) OU processes, [8] whose asymptotic 
behavior directly follows from the Fourier transform of p(x,t). 

In contrast to the standard Gaussian case (folk theorem: all solutions of a given 
Fokker-Planck equation have a common asymptotic behavior), there are no general 
mathematically rigorous statements about the asymptotic behavior of solutions of 
a general fractional Fokker-Planck equation. Specifically, this comment pertains to 
whether (if at all) and how (detailed convergence estimates) invariant (stationary, equi- 
librium) densities are asymptotically approached in the course of a nonlinear (symmet- 
ric) stable process. Anyway, for monomial drifts, analytic forms of associated invariant 
densities were explicitly elaborated in the presence of Cauchy noise in Refs. pm CGI [13] . 

It is also well known that, given the drift b potential V(x), (1), an invariant density 
for a confined Levy flight does not show any connection with the Gibbsian exponential 
of the form (3). Hence, the usefulness of the transformation (4) comes under scrutiny 
in the presence of Levy noise. 

3.2.2 Feynman-Kac (topological) route 

It is a possible to account for external perturbations by means of Levy-Schrodinger 
Hamiltonians, [15] and [5], [6] where a potential function V(x) appears as a necessary 
ingredient. Assuming that its functional form guarantees that 

= A|A|^ 2 + V. (25) 

is self-adjoint and positive in a suitable Hilbert space, we can consistently introduce the 
Levy-Schrodinger semigroup exp(— Iff^) and the fractional relative of the generalized 
diffusion equation: 

d t 9* = -\\A\»' 2 e,-V9*. (26) 
The time-adjoint equation (compare e.g. (15)) has the form 

d t B = X\A\" /2 e + Ve. (27) 

Clearly, we have here reproduced the general theoretical framework of the Schrodinger 
boundary data problem, where 9*(x,t)9(x,t) = p(x,t) stands for a probability density 
of an associated Markov process. 

Let p*(x) be a stationary (invariant, equilibrium) probability density of the perti- 
nent process. With 9(x,t) = 9(x) = exp[$(a;)], we can mimic the trial decomposition 
ansatz, Eq. (4) 9* = p exp(— $). 
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If we set exp[$(x)] = pj (x), then Eq. (27) takes the form of the compatibility 
condition, akin to that of Eq. (14): 

I AIM/2 1/2 

V = -X 1 -^-. (28) 
p* 

This identity should be compared with Eq. (8) in Ref. [17] . where an analogous 
effective potential (up to a systematic sign difference) was deduced for the fractional 
Levy-Schrodinger type equation, in the study of Levy flights in inhomogeneous media. 

In view of (26) and (27) we have a continuity equation with an explicit fractional 
input 

d tP = 9d t 9* = -A(exp $)|A|^ 2 [exp(-$)p] - V • p = -Vj . (29) 

Let us make cosmetic changes $ — > —V/2kBT, c.f. Eqs. (3) and (17) for comparison). 
Next we take advantage of p = 9*9 with 9 = exp(— (3V/2) and Eq. (28). By setting 
A = 1 and (3 = 1/ksT, we give Eq. (35) a familiar form of the transport equation 
previously introduced in a number of papers: 

d t p = -exp(-/3V/2) |A| /i/2 exp(/5V/2)p + pexp(/3\//2)|A|' l/2 exp(-/?V/2) , (30) 

c.f. formula (6) in [17] . formula (5) in [18] and formula (36) in [TJ]. There, the 
investigated process was named a topologically induced super- diffusion. We point 
out a systematic sign difference between our |A|^/ 2 and the corresponding fractional 
derivative A^ 1 / 2 of [T4l[T7] [T8]. Graphically these symbols look similar, but have different 
origin. 

The major observation at this point is, that topological Levy (specifically, Cauchy) 
processes have been embedded into a well developed mathematical framework of Ref. 
[6]. Therefore, we can reconcile any specific choice of 9 ~ p^J 2 with minimal require- 
ments upon the properties of V, Eq. (28), that guarantees the existence of Levy- 
Schrodinger semigroup and thence of the fractional transport equations (29), (30). 

3.2.3 A discord and its analysis 

The puzzling point, raised in [TJ], is that for Levy processes in external force fields, 
the Langevin approach yields a continuity (e.g. fractional Fokker- Planck) equation in 
a form 

dtP = -v(-^p) -X\A\^ 2 p (31) 



m(3 

that is very different from the previous fractional transport equation (29). 

The conclusion of Refs. [H]-[T8] was that, while assuming $ ~ —V where V is the 
external force potential (up to inessential factors), the two transport equations (29) 
and (31) are plainly incompatible. Eq. (31) seems not to correspond to any standard 
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Langevin equation with Levy noise term and b = — W/m/3 as its deterministic part 
and in reverse. 

Apart from the verbal statement this puzzling discrepancy has not been explored 
in more depth. We shall partly fill this gap in below. 

Let us take the fractional Fokker-Planck equation (31) for granted. Assume that 
p(x,t) — > p*(x), as t —>■ oo. Let the invariant probability density of the fractional 

1 /2 

equation (31) determines $ through exp[$(x)] = pj (x). We know that in case of 
Levy noise p*(x) does not have a Gibbsian structure p* ~ exp(— V/ksT), with V being 
the drift potential. 

The general problem to be addressed is: 

(i) choose a functional form of V(x) and thus the drift of the Langevin-type process, 

(ii) infer an invariant density p* that is compatible with the fractional Fokker-Planck 
equation (31), 

(iii) given p*, deduce the corresponding Feynman-Kac (e.g. dynamical semigroup) po- 
tential V by means of (28): the two dynamical scenarios (29) and (31) would thus share 
a common stationary density, 

(iv) use V in (28), (29) and verify whether the "topologically induced dynamics" is 
afline, if at all, to that associated with (31) (and thus to the underlying Langevin 
equation with Levy noise), 

(v) check an asymptotic behavior of p(x, t) in both scenarios (29) and (31) to find pos- 
sible differences in the speed (convergence time rate) with which the common invariant 
density p*(x) is approached. 

(vi) repeat the procedure in reverse by starting from (iii) and then deduce the drift 
for the Langevin equation with Levy noise; next compare the dynamical scenarios (29) 
and (31) for any common initial probability density. 

In below, we shall mostly concentrate on the above points (i)-(iii). However, their 
reverse (vi) will receive some attention as well. We note that problems (iv)-(vi) need 
more detailed computer-assisted analysis, that is postponed to a future publication. 

4 Processes induced by Cauchy noise: Invariant 
density vs semigroup (Feynman-Kac) potential 

In view of serious technical difficulties, we shall not attempt a fully fledged solution 
to the just formulated problem, for any symmetric stable jump-type process and any 
conceivable drift. Instead, we turn our attention to situations where explicit functional 
forms of invariant densities are available. 

Most of them were inferred in connection with Cauchy noise, [HIT], [8]- [13]. In par- 
ticular, attention has been paid to confining properties of various monomial drifts upon 
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the Cauchy noise. On the other hand, Levy flights through a "potential landscape" 
(topological processes of [2]- [IB]) were interpreted as (enhanced) super-diffusions and 
an issue of possible asymptotic densities has not been significantly developed. 

For a pseudo-differential operator |A| M / 2 , the action on a function from its domain 
is greatly simplified, in view of the properties of the Levy measure v^dx). We have 

(|Ar /2 /)W = - / [f(x + y) - f{x) - V -^4-\ v»(dy) 

(lAl^m*) = - I [/(* + V) - f(x)K(dy) . (32) 



The Cauchy-Levy measure, associated with the Cauchy semigroup generator |A| 1//2 = 
|V|, reads 

v\/2{dy) = (33) 
7T y 2 

By changing an integration variable y — > z = x + y, we give Eq. (38) the familiar 
form 

(|V|/)(x) =-i / fi f - f W dz (34) 

7T J \Z — Xf 

where 1/it\z — x\ 2 has an interpretation of an intensity with which jumps of the size 
\z — x\ occur. Once we set / = exp$, the "topologically induced" jump-type process 
of Refs. 1 141- 1 181 is in fact obtained. 



4.1 Ornstein-Uhlenbeck-Cauchy process 

In case of the Ornstein-Uhlenbeck-Cauchy (OUC) process, the drift is given by b(x) = 
— jx, and an asymptotic invariant density associated with 

d tP = -A|V|p + V[( 7 x)p] (35) 

reads: 

7T(cH + X 2 ) 

where a = A/7, c.f. Eq. (9) in Ref. [?]. 

Note that a characteristic function of this density reads —F(p) = — a\p\ and gives 
account of a non-thermal fluctuation-dissipation balance. The modified noise intensity 
parameter a is a ratio of an intensity parameter A of the free Cauchy noise and of the 
friction coefficient 7. 

For the Cauchy random variable X t we have E[exp(ipX t )} = exp(tA|p|) where A 
stands for an intensity of the Cauchy process. The corresponding (time-dependent) 
probability density has the form (36) with a ~ tX, e. g. p(x,t) = Xt/7i[(\t) 2 + x 2 }. 
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Here, a and likewise tX play the role of scale parameters which specify the half- 
width of the Cauchy density at its half-maximum. Since tX grows monotonically, the 
corresponding free Cauchy noise probability density flattens and its maximum drops 
down with the flow of time. 

In view of a = A/7, the frictional drift — 7X may stop the "flattening" of the prob- 
ability distribution and stabilize the density at quite arbitrary shape (with respect to 
its maximum and half- maximum related half-width), by manipulating 7. For example, 
7 > 1 would induce a significant shrinking of the distribution p*, if compared to the 
reference (free noise) probability density at any time t ~ 1/A. In parallel, a maximum 
value of the density would increase: 1/nX — > j/nX. 

Clearly, large friction has a confining effect on Cauchy noise. Confined Levy flights, 
and specifically confined Cauchy flights, were analyzed before in [10]- [13] with the aim 
to produce explicit invariant probability densities in the presence of external (confining) 
forces. Their properties proved to be quite interesting for monomial and polynomial 
choices of V(x). 

To deduce the Feynman-Kac potential V for the OUC process, given an invariant 
density p*, Eq. (36), we need to evaluate 



V{x) 

In the notation a 
reads, (2T 

A 

7T 



A 



7T 



a 2 + x 2 f/ 2 



a 2 + x 2 , b 



1 



■s/a 2 + (x + y) 2 Vcr 2 + x 2 



dy 

y 2 



(37) 



2x, R{y) = o 2 + [x + y) 2 the indefinite integral 



y 2 ^W) y 2 



dy 



X 

71 



yyja 2a \ 2a\y\ 



(3* 



Because of the singularity at y = 0, we must handle the integral in terms of its principal 
value, i.e. by resorting to f — > f_ e + / e + °°, and next performing the e — > limit. 
Taking into account that Arshx = ln(x + y/l + x 2 ), [26J, we ultimately get 



V(x) = - 

71 



2 x, 
~P + - In 

'a a 



a + x 



a — x 



(39) 



V(x) is bounded both from below and above, with the asymptotics (2/|x|)ln|x| at 
infinities, well fitting to the general mathematical construction of (topological) Cauchy 
processes in external potentials, [6]. 

Accordingly, we know for sure that there exists a topological Cauchy process with 
the Feyman-Kac potential V, Eq. (39), whose invariant density coincides with that for 
the Langevin-supported OUC process. 

At the moment, we have nothing to say about a detailed time- dependent behavior of 
the topological process and a particular scenario of an approach towards the invariant 



13 



density (equilibrium) in the large time regime. The two considered jump- type pro- 
cesses, whose dynamics is embodied respectively in the fractional Fokker-Planck equa- 
tion and in the Levy-Schrodinger semigroup (topological case) dynamics definitely stay 
in affinity, since they share a common invariant density. In the near-equilibrium regime, 
any dynamical distinction between these motion scenarios, becomes immaterial. How- 
ever, their detailed dynamical behavior far-from-equilibrium might be different and 
this issue deserves further exploration. 

There is no faithful Langevin-type representation of a topological process and in re- 
verse, even though an invariant density is common for both. Nonetheless, we have 
demonstrated that by staring from a common initial probability density, the two 
(Langevin and dynamical semigroup) motion scenarios may, in principle, end up at 
a common invariant density. 

4.2 Confined "topological" Cauchy process 

Neither the OUC process nor its topological relative are confined. For the Cauchy 
density, the second moment is nonexistent. We shall verify the outcome of the OUC 
discussion for Cauchy-type processes whose invariant densities admit the second mo- 
ment. Let us consider the quadratic Cauchy density: 

"• [x) = ljrh?- (40) 

The action of |V| upon this density can be evaluated by recourse to the free Cauchy 
evolution. 

We note that (l/y/27r)pV 2 = + x 2 ) actually is the Cauchy probability 

1 /2 

density. Let us consider f{x) = pj as the initial data for the free Cauchy evolution 
d t f = A|V|/. This takes f(x) into 

w^flwW' < 41 > 

Since 

\\V\f = -limd t f (42) 

we end up with 

]im tl0 d t f x 2 -I 

V(x) = (x) = . (43) 

The shape of this potential is quite inspiring. A minimum —A is achieved at x — 0, 
V = occurs for x — ±1, a maximum +A is reached at x — > ±oo. 

The potential is bounded both from below and above, hence can trivially be made 
non-negative (add A). Therefore, the invariant density (40) is fully compatible with the 
general construction of the Cauchy-Schrodinger semigroup and the induced jump-type 
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process, c.f. Corollary 2, pp. 1071 in [5J. This topological Cauchy process is induced 
by the Cauchy generator plus a potential function V given by Eq. (43). c.f. Corollary 
2, pp. 1071 in [6j. The process is of the jump-type and can be obtained as an e j, 
limit of a step process, e.g. jump process whose jump size is bounded from below by 
e > but unbounded from above. 

In connection with the problem (vi) o Section 3.2.3 let us note that if the quadratic 
Cauchy density (40) would actually stand for a stationary density of the fractional 
Fokker-Planck equation with a drift Eq. (31), then we should have: 

%>. = = -V(&p„)-7|V|p„. (44) 
Therefore the drift function, if any, may be deduced by means of an indefinite integral: 

b( x ) = ^— [ (\V\p*){x) dx . (45) 



If we equate 0, otherwise arbitrary integration constant, we have associated with 
p*(x) = 2/7r(l + x 2 ) 2 an admissible drift function: 

b(x) = -- (x 3 + 3x) . (46) 
8 

Thus, if we wish to deal with the Langevin process associated with the quadratic 
Cauchy density (40), the proper drift form is given above. 

4.3 "Topological" Cauchy family 

We may consider various probability densities as trial ones. Let us pay attention to 
a broader class of densities that bear close affinity with the Cauchy noise. With a 
given continuous probability distribution p we associate its Shannon entropy S(p) = 
— f p hipdx, [25]. If an expectation value _E[ln(l + x 2 )} is prescribed (e.g. fixed), the 
maximum entropy probability function belongs to a one-parameter family 

. . _ r(a) 1 
P * {X) = ^T{a - 1/2)) {l + x 2 Y [ V 

where a > 1/2, [25]. 

The gamma function r(a) = J °° exp(— t) dt we specialize to integer a = n + 1- 
values, with n > 0. ThenT(n+l) = n! and V{a-l/2) -»• r(ra+l/2) = [(2n)\^]/n\2 2n . 

The Cauchy distribution is a special case of the above p* that corresponds to a = 1. 
The density (40) is the second, a = 2, member of the a-integer hierarchy (we presume 

(7 = 1). 

To elucidate an intimate connection with the underlying Cauchy noise, let us invoke 
the Fourier transform J^fp*], e.g. a characteristic function of p*. Namely, we have, [26] : 



A +7 )l ~ {n-l)\ <9z«-i 



(4J 
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For the (unnormalized, a = 1) Cauchy density l/(a; 2 + 7 2 ), we infer (ir/j) exp(— 7|p|). 
For the quadratic (a = 2) case, we obtain (7r/27 2 )[|p| + (I/7)] exp(— "f\p\). 

Given an invariant probability density of the form (47). Assume that V(x) can be 
inferred, e.g. exists, and additionally fits the restrictions of Corollary 2 in Ref. [6]. 

Then, we can be sure that the corresponding "topological" Cauchy-type process 
can be consistently defined. For each concrete V(x), the resulting Markov stochastic 
process of the jump-type, that is determined by the Cauchy generator plus a suitable 
potential function, appears to be unique. Specifically, let us consider 

07r (1 + ar) 4 

By evaluating principal value integrals, we end up with the following expression for the 
Feynman-Kac (semigroup) potential: 

V ( x ) = onl ^ ^ + 6x2 ~ 3) • (50) 
2(1 + x z ) 

The potential is bounded from below, its minimum at x = equals —37/2. For large 
values of \x\, the potential behaves as ~ ( r y/2)x 2 i.e. shows up a harmonic behavior. 

Apart from the unbounded- ness of V(x) from above, this potential obeys the min- 
imal requirements of Corollary 2 in Ref. [6]: can be made positive (add a suitable 
constant), is locally bounded (e.g. is bounded on each compact set) and is measurable 
(e.g. can be arbitrarily well approximated by means of sequences of step functions). 

We can readily address the problem (vi) o Section 3.2.3. Be employing the density 
(49), we get, [26]: 

b( x ) = -— (5x 6 + 21a; 4 + 35x 2 + 35) . (51) 
16 

This a bit discouraging expression shows a linear friction b ~ — x for small x and a 
strong taming behavior b ~ — x 7 for large x. 



5 Conclusions and prospects 

Insightful, explicitly solvable models are scarce in theoretical studies of Levy flights, 
in the presence of external potentials and/or external conservative forces. Therefore, 
our major task was to find novel examples, that would shed some light on apparent 
discrepancies between dynamical patterns of behavior associated with two different 
fractional transport equations (29) and (31) that are met in the literature. 

We have demonstrated that so-called topological Levy processes form a subclass 
of solutions to the Schrodinger boundary data problem. The pertinent dynamical 
behavior stems form a suitable Levy-Schrodinger semigroup. The crucial role of the 
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involved Feynman-Kac potential has been identified. We have deduced these potential 
functions explicitly in a number of cases. 

The major gain of those observations is that a mathematical theory of Ref. [6j tells 
one what are the necessary functional properties of admissible Feynman-Kac potentials. 
Their proper choice makes a topological Levy process a well behaved mathematical 
construct, with a well defined Markovian dynamics and stationary density. That gives 
clear indications towards any explicit computer-assisted modeling. 

Our focus was upon confinement mechanisms that tame Levy flights to the extent 
that second moments of their probability densities are admissible. We have proved 
that the pertinent dynamical patterns of behavior stay in close affinity in the near- 
equilibrium regime and, in each considered case, admit common for both stationary 
density. In turn, this density determines a functional form of the above mentioned, 
dynamical semigroup-defining, potential function. 

Our statement of problems involved can be found in Section 3.2.3. We have analyzed 
problems (i)-(iii) and (vi). A number of important and interesting topics has been 
left aside in the present paper. Namely, one may choose a priori a simple, physics 
motivated, functional form of the potential V (like e.g. V ~ x 2 or ~ x 4 ) and thence the 
related Levy and specifically Cauchy semigroup. The major technical problem (to be 
circumvented by a computer-assisted exploration of the issue) is to deduce an explicit 
functional form of the corresponding asymptotic invariant density and to visualize 
dynamical patterns of behavior in the vicinity a stationary state. This topic has not 
received attention in papers on topological super- diffusions [L4"]-[TT]. 

Since, the Schrodinger boundary data problem allows for a construction of an inter- 
polating Markovian processes between any two a priori prescribed probability densities, 
it is of interest to fix an initial density and choose an invariant density as an asymptotic 
(terminal) datum. A comparison of a detailed temporal behavior of the Langevin-based 
and affine topological process (both sharing the chosen invariant density) is of interest. 
This issue is postponed to the future publication. 
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